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Cover:  An  airborne  radio-echo  sounding  record,  obtained  using  a 60-MHz 

sounder  from  the  ice  sheet  in  the  vicinity  of  Byrd  Station,  Antarctica 
(80° S,  120'W).  S and  b represent  the  surface  and  bottom  of  the  ice 
sheet,  and  the  horizontal  white  lines  between  them  are  radar  returns 
from  internal  layers  within  the  ice  sheet.  The  ice  thickness  of  2164 
m measured  by  the  core  taken  at  this  location  is  in  good  agreement 
with  the  radar  computation  using  the  travel /time  scale  below  the 
photo  and  a velocity  of  radiowaves  in  ice  of  168  m/psec  (two-way 
travel  time).  (This  photo  was  provided  by  Dr.  D.J.  Drewry 
of  the  Scott  Polar  Research  Institute,  Cambridge,  England.) 
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The  use  of  radio-echo  sounding  records  to  indicate  the  presence  of  internal  layers  within  large  ice  sheets  is  of  interest 
to  glaciologists  because  it  offers  a means  of  tracking  the  internal  properties  of  the  ice  sheets  over  large  distances.  The 
interpretation  of  the  reflections  obtained  in  this  manner  is  more  valuable,  however,  if  a physical  property  change  re- 
lating to  the  glaciological  regime  can  be  related  to  the  dielectric  property  change  producing  the  radio-echo  reflections. 
In  this  report,  we  use  the  measured  physical  properties  of  core  to  bedrock  taken  at  Cape  Folger,  East  Antarctica 
C[66°22'S^1 1 flE,  324-m  depth),  to  compute  a profile  of  dielectric  properties  and  from  this,  a depth-reflection 
coefficient  profile  for  comparison  with  the  observed  radio-echo  reflections.  The  measurements  available  on  physical 
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20.  Abstract  (cont'd) 

^properties  are:  density  variations,  bubble  size  and  shape  changes,  and  crystal  fabric  variations.  The  depths  of  the 
strong  reflections  shown  on  the  available  radio-echo  records  are  in  reasonable  agreement  with  the  depths  correspond- 
ing to  the  highest  reflection  coefficients  computed  from  the  combined  physical  property  measurements.  In  calcula- 
tions to  differentiate  the  separate  effects  of  different  physical  properties,  it  appears  that  density  variations  account 
for  the  primary  contributions  to  the  calculated  dielectric  property  changes  corresponding  to  the  highest  reflection 
coefficients.  However,  bubble  changes  alone  can  also  account  for  reasonable,  though  lower,  reflection  coefficients 
at  the  appropriate  depths.  Crystal  fabric  variations  correspond  poorly  with  the  reflection  locations.  Density  varia- 
tions are  normally  associated  with  depositional  events  in  the  history  of  the  ice  sheet.  However,  the  close  correspond- 
ence between  the  depths  of  the  bubble  shape  changes  (which  are  definitely  deformational  features),  and  the  depths 
of  the  density  variations,  and  between  both  of  these  and  the  radio-echo  layers,  indicates  that  deformational  events 
in  the  ice  sheet’s  history  are  represented  by  the  variations  in  the  physical  property  and  associated  radio-echo  records. 


Unclassified 


SCCuniTY  CLASSIFICATION  OF  THIS  PM€(W»nn  DM* 


PREFACE 


This  report  was  prepared  by  T.E.  Keliher,  Research  Associate,  Department  of  Meteorology, 
University  of  Melbourne,  Parkville,  Victoria,  Australia,  and  S.F.  Ackley,  Research  Physicist,  Snow 
and  Ice  Branch,  Research  Division,  U.S.  Army  Cold  Regions  Research  and  Engineering  Laboratory 
(CRREL). 

The  study  was  conducted  under  DA  Project  4A161 102AT24,  Research  in  Ice  and  Frozen  Ground; 
Task  A3,  Research  in  Terrain  and  Climatic  Constraints,  E 1 , Cold  Environment  Factors;  Work  Unit 
002,  Adhesion  and  Physics  of  Ice. 

Dr.  Steven  A.  Arcone  and  Dr.  Anthony  ).  Gow  of  CRREL  reviewed  the  technical  content  of 
this  report. 

The  authors  thank  Dr.  W.F.  Budd,  Antarctic  Division,  Department  of  Science,  Australia,  for  pro- 
viding the  unpublished  data  on  ice  core  physical  properties  used  in  this  report  and  for  encouraging 
and  supporting  this  work.  Both  of  the  authors  also  gratefully  acknowledge  the  support  of  the 
Meteorology  Department,  University  of  Melbourne,  Dr.  U.  Radok,  Chairman,  during  this  investi- 
gation. S.F.  Ackley  particularly  acknowledges  the  support  of  CRREL. 


ACCESSION  lor 

NTIS  White  Section  & 

PDC  Boll  Section  □ 

UfMNNOUNCn)  □ 

iUSTI'ICATION  . . 


DISTRIBUTION/ AYA2-. WHllTY  COPES ] 

0."  SP.CIAL 


I 


r 


'-'T 


n 


CONTENTS 

Page 


Abstract * 

Preface Hi 

Nomenclature v 

Introduction 1 

Derivation  of  the  power  reflection  coefricient-depth  variation t 

Possible  sources  of  dielectric  constant  variation  with  depth  in  ice  sheets  2 

Ice  density  variations 3 

Effect  of  macroscopic  impurity  layers 3 

Effect  of  microscopic  impurity  layers 3 

Changes  in  inclusion  geometry 4 

Effect  of  crystal  orientation  changes 5 

Physical  properties  of  the  Cape  Folger  core 6 

Density-depth  profile 7 

Air  bubble-geometry  depth  profile 7 

Crystal  orientation-depth  profile  8 

Radio-echo  sounding  data 8 

Results  and  conclusions 9 

Literature  cited  12 


ILLUSTRATIONS 


Figure 

1 . Reflection  amplitude  as  a function  of  layer  thickness  measured  in  wavelength  of  the 

incident  wave  2 

2.  Fabric  diagram  from  171-m  depth  at  Cape  Folger 6 

3.  Distribution  of  c-axes  from  the  vertical  at  171-m  depth  obtained  from  the  fabric  diagram 

in  Figure  2 averaged  over  3°  intervals 6 

4.  Variation  of  ice  density  with  depth  at  Cape  Folger  and  the  results  of  a linear  regression 

of  density  upon  depth 7 

5.  Variation  with  depth  in  the  ratio  of  long  axis  to  short  axis  a/6  of  air  bubbles  from  the 

Cape  Folger  core  7 

tt/2 

6a.  First  moment  of  the  distribution  ff  = V 0n(0) &0  of  c-axcs  from  the  vertical 8 

6b.  Variation  of  the  dielectric  constant  caused  by  mean  crystal  fabric  angle  variations 

shown  in  Figure  6a 

7.  Dielectric  constant  and  power  reflection  coefficient  vs  depth 

8a.  The  quantity  elb  and  the  corresponding  reflection  coefficient  vs  depth  with  and  with- 


out crystal  effects 10 

8b.  The  quantity  eb  and  corresponding  power  reflection  coefficient  vs  depth  as  in  Figure  8a  10 
9a.  The  dielectric  constant  and  reflection  coefficient  for  the  various  cases  of  bubble  con- 
figurations as  a function  of  depth  using  the  linear  regression  density  1 1 

9b.  The  dielectric  constant  and  reflection  coefficient  for  crystal  fabric  effects  as  a function 

of  depth  using  the  linear  regression  density 1 1 


\ 


iv 


vO  00 


NOMtNlLATURt 


*f|  depolarization  factor 

1 1,  h,  i diameters  in  v,  y,  i space  of  a general  ellipsoid 

d,  thickness  of  laser  i (metets) 

f radio  frequency  (ft/) 

/ v^T 

t thickness  of  an  ash  or  dust  laser 

n total  number  of  layers 

P proportion  of  ash  to  ice  in  a layei 

r,  tM  reflection  amplitude  ol  layer  i 

K powei  reflection  coefficient  ot  mixture  ol  ash  and  ice 

Kt  powei  reflection  coefficient  ot  laser  / 

s vai  iable  lelating  general  ellipsoid  to  the  "standard"  ellipsoid 

Vj  volume  density  ot  component  2 

bulk  impedance  ot  layer  i 

sur  lace  impedance  equivalent  to  hall  space  impedance  below  layer  i 
>(  propagation  constant  of  layer  / 

e()  permittivity  ol  tree  space 

et  relapse  dielectric  constant  (to  tree  space)  of  laser  i 

e’(co)  dielectric  constant  ot  a medium  as  a function  of  frequency 

«?i . r_.  dielectric  constants  ot  the  components  1 and  2 ot  a mixture  dielectric 

t total  dielectric  constant  ot  the  mixture 

tt  dielectric  constant  ot  a mixture  containing  ellipsoids  with  /.  field  paiallel  to  the  "long" 
axis  it 

tp  dielectric  constant  of  a mixture  containing  ellipsoids  with  /'•field  paiallel  to  the  "short" 
axis  t> 

t jp  dielectric  constant  of  a mixture  containing  ellipsoids  with  a random  spatial  orientation 

(%  dielectric  constant  ot  a mixture  containing  spheres  as  the  included  second  dielectric 

r,  dielectric  constant  at  infinite  frequency 

Aen  contribution  to  the  dielectric  constant  from  polarization  process n 

dielectric  constant  of  ice  measured  perpendicular  to  the  c axis 
1 1)  dielectric  constant  of  ice  measured  parallel  to  the  c -axis 
f’  real  pait  ot  the  complex  dielesttic  constant 

e"  imaginary  pait  of  the  complex  dielectric  constant 

0 angle  measured  in  the  veitical  between  l-M  propagation  direction  and  crystal  axis 

A wavelength  of  the  incident  radiation 

9 mean  vertical  angle  ot  thec-axis  distribution 


i 

I 

<p  azimuth  ingle  of  crystal  axis 

pL  (/)  linear  regression  of  density  upon  depth 
p,<f  density  of  ice 

p,  density  of  layer/ 

Hq  magnetic  permeability  of  free  space  (henrys/m) 

0,  conductivity  of  layer  i 

r„  relaxation  time  of  the  polarization  process  n 

co„  angular  frequency  (n)  of  the  polarization  process/? 

u angular  frequency  = 2wf  (sec'' ) 

tan  6 loss  tangent  of  a medium  = e“/e' 

dll  sin  OdtidQ 

percenuge  of  ice  crystals  with  crystal  axes  at  angles  6,  p 


* COMPARISON  BETWEEN  DERIVED  INURNAL  DIELECTRIC 
PROPERTIES  AND  RADIO-ECHO  SOUNDINC.  RECORDS 
OP  TNI  ICE  SHEET  AT  CAPE  FOLC.ER,  ANTARCTICA 

I I kelihei  and  Yl  Ackley 


INTRODUCTION 

1 ho  use  of  radio-echo  sounding  records  to  indicate 
the  piesence  of  internal  lasers  vc ilbm  large  ice  sheets 
is  of  interest  to  glaciologists  because  it  oftet'  a means 
of  Hacking  the  internal  properties  of  the  ice  sheets 
osei  large  distances.  The  glaciological  interpretation  is 
difficult,  hoscever,  because  the  cause  of  the  internal 
laxeitng  is  not  cleat  In  tact,  seseial  physical  causes  loi 
the  dielectnc  changes  of  the  ice  producing  internal  te- 
tlections  have  been  otteiesf  (Harrison  Id7.1|  lach  of 
these  causes  is  glatiologically  significant,  but  in  a dil 
ferent  mannei.  so  that  large  scale  internal  lasermg 
within  the  ice  cannot  be  unequivocally  mteipieted  as 
a partnulai  stratigraphic  horizon  relating  to  a chemi 
sal.  detormational,  or  temper  atuiedepositional 
characteristic  of  the  ice.  It  is  also  possible  that  a sin- 
gle ladai  reflection  mas  arise  from  the  sector  sum  of 
reflections  from  seseral  lasers. 

In  this  irpotl  we  first  review  a computational 
method  to  obtain  reflection  coefficients  (l  inlor  and 
luacek  1 A)  and  next  testess  seseia!  ps'ssible  phssical 
mechanisms  that  san  pioduce  the  necessary  dielectric 
constant  variations  leading  to  reflections,  he  then  use 
the  measuied  phssical  pioperties  of  th-  core  to  bed 
rock  tromC'ape  folgei  neai  Cases  Station,  fast 
Antarctica  (bb'’.?0'S.  1 1 l*E;  324-m  depth),  to  calcu- 
late the  contribution  ts>  the  dielectric  properts  change 
with  depth  trssm  each  s'bserved  phssical  cause,  f rout 
the  calculated  dielec  Hie  profile,  we  compute  a derived 
power  letlectnsn  coefficient-depth  profile.  We  then 
ci'mpore  this  computed  letlevtis'n  coefficient  profile, 
s'btamed  trs»m  the  individual  and  combined  phssical 
properts  changes  of  the  ice,  with  the  actual  radio- 
echo  sounding  reflection  profile  in  an  effort  to  inter 
pret  the  causes  of  the  internal  reflections. 


DERIVATION  OF  THE  POWER  REFLECTION 
COEFFICIENT-DEPTH  VARIATION 

Ihe  powei  reflection  coefficient  is  calculated  in  the 
following  mannei  attei  Lmlot  and  litacek  (l*>75)  and 
hail  (1958).  Ihe  characteristic  bulk  impedance  /,  of 
las  er  r is  defined  bs 

(,) 

ri 

wheie  ; T v - ' 

ftp  ” magnetic  peimeabilits  of  free  space 
co®  angular  ftequencs  'ft  (/  in  Hz) 

T,  * piopagation  constant. 

Generally  y,  is  a complex  quantity  defined  bs 

7?  * e0  Ro  *i w*  + tf*o  4,i w (*) 

where  e,j  • peimittivits  of  free  space 

f,  « relative  lossless  dielectric  constant  of  the 
medium 

ot  ■ conductivity. 

In  our  formulation  the  ice  medium  is  consideied  loss 
less,  sc'  we  only  conxidei  the  leal  pait  ot  the  piopagation 
constant,  or 

r,  co\70n0^  . (T) 

Paien  and  Robin  (Id'S)  have  considered  the  effects  of 
conductisits  changes  in  ice  sheets  (corresponding  to  a 
change  in  the  imaginais  pait  of  eq  2)  on  internal  reflec- 
tions but  have  indicated  that  these  effects  aie  signifi- 
cant at  depths  over  1000  m deeper  than  the  core 


r 


N 


considered  here.  We  discuss  this  aspect  further  in  a 
later  section 

The  reflection  amplitude  r(  it  a boundary  wilhin 
the  medium  it  given  by 


(4*) 


and  the  reflection  power  (wh.it  iv  detested  hv  the  re- 
ceiver) it 


i 


r: 


(4b) 


/,_!  iv  the  characteristic  bulk  impedance  ot  the 
/-Ivt  lav er  as  defined  above,  while  iy  iv  the  surface 
impedance  of  the  /th  laser  and,  after  Wait  1 1 **58).  ' 
defined  bv 
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£ K\  */, tjnh  >,<*, 

' tjnh^d, 


(5) 


Here,  d,  ■ the  thickness  of  layer  / 4nd  y,  ■ the 
propagation  constant  4v  given  above  1 he  bottom  or 
nth  laver  at  bedrock  iv  taken  to  be  vemi-infinile  and 
therefore 
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I he  procedure  to  obtain  the  estimated  reflection 
coefficient  profile  with  depth  iv  to  "build  up"  the  ice 
sheet  from  the  bottom,  laser  b\  laser,  using  some  esti- 
mate ot  the  propagation  constant  based  on  the  ice 
properties  at  each  depth.  The  bulk  impedance  of  a 
laser  is  calculated  from  ecj  I . tq  5 is  then  used  with 
the  calculated  bulk  impedance,  estimated  propagation 
constant  and  estimated  laser  thickness  to  obtain  the 
surface  or  input  impedance  of  the  top  which  replaces 
the  "lasered"  dielectric  below  by  a homogeneous  half 
space  with  equivalent  surface  impedance. 

The  variation  of  reflection  coefficient  with  depth, 
which  is  the  source  ot  internal  reflections  within  the 
ice  sheet,  iv  therefore  dependent  on  the  "mismatch" 
m impedance  between  the  bulk  impedance  of  the  laver 
under  studs  and  the  surface  impedance  of  the  half  space 
below  the  laser  From  eq  I and  5 we  see  that  these  two 
quantifies  are  dependent  onls  on  the  dielectric  constant 
and  thickness  ot  the  lasers  for  a given  tadar  wavelength. 


POSSIBLE  SOURCES  OF  DIELECTRIC  CONSTANT 
VARIATION  WITH  DEPTH  IN  ICE  SHEETS 


In  this  analysis  one  effect  we  do  not  consider  is  the 
contribution  of  interference  between  reflected  waves 


— '■* 


f iguir  /.  Ketlec tkm  amplitude  as  a 
function  of  layer  thickness  measured 
in  sssiselength  of  the  incident  wuve 
(after  link*  end  lirocek  1975).  ,4s 
shows?  here,  the  curses  represent  an  ice 
sheet  floating  on  *»  atei.  Similar  curses 
can  he  generated  tor  analogous  3-layer 
media  /e.g.,  air-snow .'he-hedrixk  (as 
discussed  here)/. 


to  the  observed  reflections  within  ice  sheets.  When  the 
thickness  of  a layei  is  either  an  even  or  an  odd  multiple 
of  a quarter  wavelength  of  the  incident  wave,  construc- 
tive or  destructive  interference  occurs  between  two 
reflected  waves.  The  reflected  power  can  then  be  en- 
hanced or  reduced  if  a laser  thickness  is  at  the  quarter 
wavelength  multiple  value.  To  accurately  evaluate  this 
effect  is  beyond  the  resolution  of  the  currently  available 
data  on  the  physical  properties  of  the  ice.  The  quartet 
wavelength  phenomenon  would  not.  a priori,  be  expected 
to  be  highly  probable  since  it  demands  an  accurate 
match  of  input  wavelength  with  the  particular  layer 
thickness  in  the  ice  sheet. 

In  the  future,  a simple  test  on  whether  a reflection 
occurs  or  disappears  because  of  a match  between  the 
wavelength  and  layei  thickness  could  be  simply  per- 
formed by  sounding  at  two  slightly  different  frequenc  ies, 
since  the  reflection  amplitude  changes  quite  rapidly 
with  a slight  mismatch  between  the  lasei  thickness  and 
the  quarter  wavelength,  as  shown  in  f iguie  1 (after 
Linlor  and  |iracek  1975).  Therefore,  if  the  quartet 
wavelength  thickness  match  were  a factor,  the  re- 
flected power  from  the  laser  under  consideration  would 
be  markedls  different  at  the  two  frequencies,  while  it 
would  not  change  significantly  if  the  quartet  wavelength 
phenomenon  were  absent.  Because  internal  reflections 
have  been  nested  since  the  advent  of  radio  echo  sounding 
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it  frequencies  varying  from  30  MH/  to  the  GHa  range, 
we  assume  initially  that  the  observed  echoes  are  de- 
rived from  depth  changes  in  ice  dielectric  properties 
(leading  to  a change  in  impedance  with  depth  as  given 
in  eqs  1, 2.  and  5)  rather  than  from  a fortuitous  coinci- 
dence of  layer  thickness  and  quarter  wavelengths. 

The  possible  mechanisms  for  dielectric  sanations  in 
ice  sheets  will  be  briefly  reviewed  here  The  mechanisms 
mentioned  have  been  (i)  ice  density  sanations  iRobin 
et  al.  1969),  (ill  macroscopic  impurity  lasers  [dust  an  t 
volcanic  ash  (Robin  et  al.  1969)  (,  (ml  microscopic 
impurity  layers  (Paren  and  Robin  1975),  (is)  a previ- 
ously unstudied  effect  for  this  application,  changes  in 
air  bubble  geometry  trom  spherical  to  ellipsoidal  shapes, 
and  (v)  crystal  orientation  changes  ; Harrison  1973) 
Except  for  microscopic  impurity  variations,  all  ot  these 
give  rise  to  small  changes  in  dielectric  permittivity . or 
the  real  part  of  the  dielectric  constant. 

Ice  density  variations 

The  ice  density  variation  is  given  as  (Gudmandsen 

1971) 

y/7  = 1 + 0.85  pj . (7) 

The  variation  in  dielectric  constant  with  density 
arises  trom  the  mixture  of  two  media,  air  and  ice,  with 
different  dielectric  constants,  (e1K  - I,eite  = 3.17). 
The  density  relationship  (eq  7)  is  empirical,  based  on  a 
fit  to  experimental  data  (Robin  et  al.  1969).  The  ice 
density  variation  and  the  effect  of  macroscopic  impuri- 
ties are  specific  examples  of  a more  general  treatment 
for  anv  geometric  dielectric  inclusion,  which  we  discuss 
later. 

Effect  of  macroscopic  impurity  layers 
(dust  and  volcanic  ash) 

The  effect  of  these  layers,  arising  from  the  fallout 
of  volcanic  eruptions,  is  similar  to  the  effect  described 
above  for  density;  namely,  the  dielectric  constant  of  a 
mixed  ice  and  ash/dust  layer  will  be  a function  of  the 
dielectric  constants  ot  ice  (3.1 7)  and  of  ash  (e.g.,  4.8) 
and  the  relative  proportions  of  the  two  components. 
After  Robin  et  al.  (1969),  the  power  reflection  coef- 
ficient of  a mixture  of  ash  and  ice  R contained  in  a 
layer  of  thickness  V is 

R = (0.4rr«/X)2  (8) 

where  P - proportion  of  ash  to  ice  ( P : 1 ) 

S » layer  thickness 
X • wavelength. 


An  example  al  the  wavelength  corresponding  to  35  MHz 
indicates  layers  I and  5 cm  thick  containing  10%  ash 
would  have  reflection  coefficients  of  - 74  dB  and  - 60 
dB,  which  are  within  the  observable  range  for  most 
current  echo  sounding  receivers.  Ash  layers  have  been 
observed  at  several  locations  within  the  deep  core  ob- 
tained at  Bvrd  Station,  Antarctica  (Gow  and  William 
son  1971),  although  the  thicknesses  are  lower  than 
those  of  the  above  example,  ot  the  order  of  0.05  cm  at 
10%  concentration. 

Effect  of  microscopic  impurity  layers 

Paren  and  Robin  (1975)  argue  that,  as  well  as  the 
changes  in  the  real  part  ot  the  dielectric  constant  given 
bv  eqs  3 and  4,  there  is  also  an  impedance  mismatch 
caused  by  variations  in  the  conductivity  or  loss  tangent, 
tan  6,  ot  the  ice  (tan  6 = e"/e’  where  e",  the  imaginary 
part  ot  the  dielectric  constant  = u/cot0).  In  their 
formulation,  the  reflection  coefficient  R at  the  inter- 
face of  a two-layer  medium  differing  only  in  loss  tan- 
gent is 

R = 1 1/4  A(tan  5)| 2 (9) 

where  the  two  layers  are  characterized  bv  loss  tangents 
of  tan  6 and  tan  5 + A (tan  5),  respectively. 

Making  use  of  the  model  of  Paren  and  Walker  (1971). 
thev  suggest  that  loss  tangent  variations  of  the  order 
necessary  to  produce  reflections  may  occur  because  of 
impurity  distributions  varying  within  the  ice,  for  ex- 
ample, impurities  either  uniformly  distributed  or  con- 
centrated at  grain  boundaries,  without  requiring  that  the 
total  concentration  vary.  This  factor  remains  fairly 
constant  with  depth,  making  it  relatively  greater  than 
the  density  variation  at  great  depths,  because  density 
fluctuations  are  smoothed  out  at  depth  bv  the  over- 
burden pressure.  Ihe  loss  tangent  effect  also  increases 
with  increasing  temperature,  which  would  enhance  its 
effect  relative  to  density  variations  in  the  deeper,  warmer 
layers  where  it  is  proposed  as  a reflection  mechanism. 

We  have  not  considered  this  mechanism  in  the  present 
analysis,  since  its  main  application,  according  to  Paren 
and  Robin  (1975),  most  likely  applies  at  a depth  greater 
than  1 500  m,  which  is  much  greater  than  the  measured 
ice  depth,  324  m,  at  Cape  Folger.  As  well,  we  note  the 
mechanism  depends  upon  a var  iation  in  taniS,  which 
has  not  been  proven  to  exist  within  the  ice  sheet  at  this 
location.  I he  meltwater  conductiv  ity  measurements 
(DC)  for  the  Cape  Folger  core  indicate  that  the  fluctua- 
tion ot  this  quantity  with  depth  is  so  small  as  to  be  ob- 
scured bv  the  measurement  error  and  possible  handling 
contamination  (V.  Morgan,  personal  communication). 


1 


3 


ll 


r I 

i 


Within  the  limits  of  these  measurements,  no  attributable 
reflection  exists  based  on  conductivity  changes  with 
depth.  Therefore,  since  density  and  other  variations  do 
exist  and  are  measurably  large  enough  to  account  for 
some  of  the  variations  in  calculated  and  observed  re- 
flection strengths,  we  do  not  feel  justified  in  invoking 
the  Paien  and  Robin  hypothesis. 

Changes  in  inclusion  geometry 

This  effect  may  be  looked  upon  as  a generalization 
of  the  density  and  macroscopic  impurity  effects  where 
the  geometry  as  well  as  the  volume  of  the  included 
dielectric  is  also  considered,  evidence  from  the  Cape 
f olger  core  indicated  substantial  deviation  in  the  air 
bubbles  from  spheroids  with  depth,  justifying  an 
investigation  into  the  possible  effects  of  these  geometric 
changes. 

Some  work  has  been  done  in  explaining  the  varia- 
tion of  snow  dielectric  properties  with  the  changing 
geometry  of  air  inclusions.  Evans  (1965)  has  sum- 
marized this  work,  which  consists  primarily  of  attempts 
to  explain  the  changes  in  snow  dielectric  properties 
using  the  empirical  "Formzahl”  concept  first  postulated 
by  Weiner  (1910).  I his  approach  is  characteristic  of 
the  earlier  work  in  heterogeneous  dielectrics,  but  be- 
cause of  its  empirical  basis,  has  certain  drawbacks. 

More  recently,  van  Beek  (1967)  has  rcpUccd  these 
empirical  formulae  by  a derivation  based  on  solution 
of  electrostatic  theory  tor  dispersed  media  of  analyti- 
cal shapes. 

I he  behavior  for  dilute  mixtures  of  varying  geome- 
tries has  been  well  verified  using  this  formulation. 

Paren  (unpublished  Ph.D.  thesis)  has  shown  that  the 
limited  data  on  snow  which  falls  somewhat  beyond  the 
"dilute"  definition  can  also  be  better  reconciled  by  a 
variation  on  this  analytical  approach  given  by  Bottcher’s 
mixing  equation.  In  that  formulation,  the  average  value 
of  the  dielectric  constant  of  the  parent  medium  is  re- 
placed by  the  total  dielectric  constant  of  the  mixture. 
This  procedure  apparently  gives  a better  estimate  of 
the  dielectric  constant  at  high  concentrations  of  the 
dispersed  medium  but  again  is  a simple  variation  on 
the  analytical  models  which  work  well  at  low  concen- 
trations. Sweeney  and  Colbcck  (1974)  have  used  the 
theoretical  formulation  in  analyzing  wet  snow  measure 
ments  and  have  found  good  agreement  between  the 
experimental  results  and  theoretic, ii  predictions  when 
the  dielectric  properties  of  the  constituents  arc  known. 

After  van  Beek  (1967),  the  effective  dielectric 
constant  of  a mixture  with  dielectric  constant  of  the 
medium  e,  containing  a dilute  mixture  of  particles  (in 
our  case,  air  or  dust)  with  dielectric  constant  e2  and 
total  volume  density  v2  i* 


e 


e,  ♦ 1-4)  (1  - v2)  ♦ir2|(f2  -c,) 
e,  ♦Aj(1  -tr2)(«2  -«|) 


(10) 


Here  v2  » (pkf  - pmMi)/picc  and  4 ( is  the  depolariza- 
tion factor  which  is  dependent  only  on  the  physical 
shape  of  the  particles,  most  generally,  ellipsoids.  After 
Stratton  (1941 ), 


A : = i-  ahe 
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ds 


2 ^ (^♦f)'/*(6**j)'/2(cJ4,)1/2(/2  4 j) 
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where  / - <z,  b,  c,  and  is  the  axis  along  which  the  £ field 
is  aligned.  It  is  also  required  that 

*a+*b**c  ■ 1 (12) 


where  a,  b,  and  r arc  the  magnitudes  of  the  axes  of  the 
ellipsoids  (x,  y,  z)  and  the  equation  of  the  ellipsoid: 


a2  + s b2  + s c2  + $ 


(13) 


where  s is  a variable  indicating  the  ellipsoid  confocal 
with  a standard  ellipsoid  whose  equation  is 


*£  +j£  +jL  a 1. 

Q2  b2  C2 


(14) 


Tor  the  cases  considered  here,  the  ellipsoids  arc  as- 
sumed to  be  prolate  spheroids,  i.c.,  b = c and  a > b. 
I he  general  expression  for  the  depolarization  factor 
then  is 


A, 


= tdzL  f ds 

2 { (j*a2)3/2(J  + 62) 


(IS) 


Along  the  minor  axis  b,  the  depolarization  factor 
from  eq  12  is. 

-4b  = J (1  ~At).  (16) 

for  the  prolate  spheroid  case,  the  solution  of  the 
integral,  given  the  ratio  a/6,  is  |from  van  Beck  ( 1 967)  | 


1 


4 


.1  . jn  u*» 


' 1 

t 


A * — — * x 

(a/6)2  - 1 [(a/6)2  - 1 1 3/2 

x «/»|(a/6M(e/6)2-1|*|.  (17) 


We  consider  four  c ases  and  now  give  the  respective 
expressions  for  the  dielectric  permittivity  for  each 
case.  These  are  obtained  by  solving  eq  10  with  the 
appropriate  value  of  A j given  by  eq  16  or  1 7. 

Case  I.  The  incident  £-field  is  aligned  along  the 
a-axes  (long  axes)  of  the  ellipsoids  (which  are  parallelly 
aligned).  Then 
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Case  2.  The  incident  £-field  is  aligned  along  the 
6-axes  (short  axes)  of  the  ellipsoids  (again  all  of  parallel 
alignment): 


undetected.  From  a theoretical  argument,  a possible 
anisotropy  in  dielectric  constant  can  be  postulated  in 
the  MHz  range. 

The  frequency  variation  of  the  dielectric  constant 
can  be  written  as  (lohari  and  Charette  197$) 


«'(«) 


* t, 


(22) 


where  em  is  the  dielectric  constant  for  frequencies 
greater  than  10,<  Hz,  and  Aen  the  contribution  from, 
and  rn  the  relaxation  time  of,  the  various  polarizations. 
ea*h  represented  by  a simple  Debye  relaxation  process. 
From  eq  22  we  see  that  if 

w2r2  >>  1 (23a) 


or 
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1 + 


e,  ♦ 4b  (e2  - «|) 


(19) 
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(23b) 


Case  3.  The  f-field  is  not  aligned  with  the  ellipsoid 
axes  (random  spatial  orientation  of  the  ellipsoids): 

* - <ab  = y «*♦  jV  (20) 


then  the  contribution  of  the  relaxation  process  Aen  to 
the  real  part  of  the  dielectric  constant  f(u>)  at  frequency 
u>  is  essentially  zero. 

The  inequality  in  eq  23  corresponds  to  the  frequency 
of  the  relaxation  cj„  = 1/rn: 


Case  4.  The  ellipsoids  are  all  spheres  (a  = 6).  For 
this  case  the  value  of  At  is  1/3  and  the  solution  for  the 
permittivity  is  the  model  solution  analogous  to  the 
empirical  equation  from  Robin  ct  at.  (1969)(eq  7). 
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For  the  Cape  Folger  case  under  study  here,  no  layer 
was  found  where  | e4b  - e,|  > 0.001 2.  The  depth  varia- 
tions between  e^,  ef  and  e4,  eb  are  larger  than  this 
value,  so  it  can  usually  be  assumed  eJb  „ e,.  Equation 
20,  instead  of  eq  21,  was  therefore  also  used  to  give 
the  dielectric  constant  for  bubbles  assumed  spherical 
as  well  as  for  randomly  oriented  ellipsoids. 


Effect  of  crystal  orientation  changes 

Recent  measurements  (Johari  and  Charette  1975) 
have  indicated  that  the  dielectric  constant  values  for 
single  and  polycrystalline  ice  agree  within  ± 0.2%  at 
measured  frequencies  of  35  and  60  MM”..  From  the 
experiments,  however,  it  was  determined  a maximum 
difference  of  1%  between  the  dielectric  constant 
perpendicular  and  parallel  to  the  c-axis  would  remain 


w„  « w.  (24) 

That  is,  if  the  frequency  of  interest  is  higher  than  the 
frequency  at  which  the  relaxation  process  occurs,  then 
the  relaxation  process  in  questiui  does  not  contribute 
to  the  dielectric  constant  at  the  higher  frequency.  For 
example,  the  well-known  molecular  orientation  processes 
that  occur  in  the  kHz  frequency  range  in  ice  do  not  con- 
tribute to  the  dielectric  constant  at  35  to  100  MHz. 

The  electronic  polarization  anisotropy  that  occurs 
in  the  infrared  region  causes  a At  between  the  c-axis 
perpendicular  and  the  c-axis  parallel  of  0.0037  at  - 10°C 
()ohari  and  Charette  1975).  Harrison  (1973)  concluded 
that  a fractional  change  of  - 10"4  in  the  dielectric 
permittivity  combined  with  multipie  layering  was  suf- 
ficient to  produce  detectable  reflections  in  radio-echo 
sounding.  The  theoretical  change  in  permittivity  of 
0.0037  caused  by  crystal  orientation  changes  could  then 
account  for  the  observed  echoes  in  the  megahertz  region 
if  sufficient  crystal  orientation  change  occurs  in  these 
layers.  This  mechanism  must  still  be  cautiously  invoked 
since  the  changes  in  permittivity  arising  from  crystal 
orientation  changes  are  very  small  and  still  have  not 
been  verified  experimentally  in  the  frequency  range  at 
which  radio-echo  sounders  operate. 
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Figure  3.  Distribution  of  c-axes  from  the  vertical  at 
1 71 -m  depth  obtained  from  the  fabric  diagram  in 
Figure  2 averaged  over  3°  intervals. 


the  vertical  angle;  i.e.,  we  assume  the  crystals  are  uni- 
formly distributed  in  a cone  over  the  vertical  angle, 
or  n(0,<t> ) = n(0)  only.  Figure  2 gives  an  example  of 
a fabric  diagram  from  Cape  Folger  indicating  this 
assumption  is  probably  not  generally  valid.  However, 
there  is  usually  more  vertical  than  azimuthal  variation 
in  crystal  fabrics  from  ice  cores;  so  the  main  features 
of  the  dielectric  constant  depth  variation  are  probably 
due  to  vertical  changes  in  orientation.  Figure  3 shows 
the  data  from  Figure  2 replotted  to  give  the  vertical 
c-axis  distribution  only. 
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The  crystal  orientation  effect  at  a given  depth  on 
the  dielectric  constant  can  be  estimated  as  follows 
(Clough  1977): 


e(/)  = + Aen-x 


/«/«  «n0  cos 4>dSl 

I,f0da 


(25) 


where  Ac,  * 0.0037,  6 is  the  angle  between  the 
propagation  direction  and  the  crystal  axis  in  the  verti- 
ral,  0 is  the  azimuthal  angle  and  dil  - sin 0d0d<t>.  For 
a first  approximation,  we  consider  only  the  changes  in 


PHYSICAL  PROPERTIES  OF  THE 
CAPE  FOLGER  CORE 

The  data  on  impurity  levels  arc  based  entirely  on 
meltwater  conductivities  and  are  judged  to  be:  (1)  at 
such  a low  level  that  we  assume  the  loss  tangent 
tan  6 < 0.01 , at  MHz  or  higher  frequencies  - i.e.,  the 
material  could  be  considered  lossless;  and  (2)  of  insuf- 
ficient variation  with  depth  to  cause  reflections  - i.e., 

A tan  6 <<  0.01.  No  volcanic  dust  or  ash  bands  were 
observed  in  the  core.  Therefore,  the  observed  reflections 
were  assumed  to  originate  from  the  dielectric  effects  of 
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Figure  4.  Variation  of  ice  density  with  depth  Figure  5.  Variation  with  depth  in  the  ratio  of  long 

at  Cape  Foiger  and  the  results  of  a linear  re-  axis  to  short  axis  a/b  of  air  bubbles  from  the  Cape 

gression  of  density  upon  depth.  Foiger  core.  Data  obtained  by  Matsuto  (unpubi) 

provided  by  Budd  ( 1972). 

the  three  mechanisms  discussed  as  follows:  (i)  density  the  points.  The  dielectric  constant  was  then  calculated 

variations,  (ii)  air  bubble  geometry  variations,  and  from  eq  20  and  the  power  reflection  coefficient  from 

(iii)  crystal  orientation  variations.  eqs  1,3,4  and  5. 

After  discussing  these  physical  property  variations 
and  the  available  radio-echo  sounding  data,  we  then  Air  bubble-geometry  depth  profile 

give  the  computed  dielectric  constant  variations  based  Figure  5,  based  on  data  measured  by  Matsuto  (un- 

on  the  individual  and  combined  physical  properties.  published),  shows  the  variation  in  bubble  geometry  with 

depth  measured  from  ice  thin  sections  taken  primarily 

Density-depth  profile  in  the  horizontal  plane;  the  figure  gives  both  the  average 

Figure  4 indicates  the  variation  of  density  with  depth.  and  maximum  values  of  a/b.  Measurements  were  taken 
The  densities  are  quite  high  even  in  the  upper  parts  of  at  the  depths  indicated.  As  mentioned  earlier,  the 

the  profile  since  the  core  was  obtained  in  an  ablation  geometry  was  assumed  to  be  prolate  spheroids  charac- 

region  of  the  ice  sheet.  The  density  profile  is  obtained  terized  by  the  major  and  minor  diameters  a and  b. 

from  5-cm  cubes  uken  at  approximately  10-m  intervals  Prolate  spheroids  are  the  conceptually  simplest  geometry 

and,  in  some  cases,  the  indicated  density  may  be  of  derived  from  the  experimental  data.  Figure  5 thus  gives 

uncertain  value.  For  example,  the  value  at  195-m  depth  the  ratio  a/b  as  a function  of  depth.  Budd  (pers.  comm.) 

appears  to  be  the  result  of  a single  measurement  and  has  indicated  that  the  portions  of  the  core  where  high 

therefore  might  be  in  error.  In  contrast,  the  relatively  a/b  ratios  are  present  (below  - 165-m  depth)  are  also 

low  density  at  225-m  depth  is  consistent  w,*h  the  characterized  by  alignment  of  the  long  bubble  axes, 

trends  immediately  above  and  below,  so,  it  is  believed  roughly  in  the  direction  of  glacier  flow  or  maximum 

it  represents  a "true”  value  for  that  point.  These  points  surface  slope.  The  dielectric  constant  was  calculated 

cause  considerable  dielectric  constant  change  and  are  of  for  the  four  cases  indicated  above  using  the  average  a/b 

particular  interest  in  the  reflection  coefficient  calcula-  taken  from  Figure  5 and  the  density  information  from 

tion  discussed  later.  For  the  purposes  of  the  dielectric  Figure  4.  Linear  interpolation  between  data  points  was 

constant  calculation,  the  density  was  obtained  from  the  made  in  both  cases  to  obtain  equally-spaced  dielectric 

data  points  indicated  or  a linear  interpolation  between  constant-depth  profiles. 
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Figure  6a.  First  moment  of 
the  distribution 
7T/2 

8 = V 0r»(0)A0 
of  c-axes  from  the  vertical. 


Crystal  Fabric  Variation 
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Figure  6b.  Variation  of  the  dielectric  constant  caused  by 
mean  crystal  fabric  angle  variations  shown  in  Figure  6a. 


Crystal  orientation-depth  profile 

Examples  of  crystal  orientation  data  from  the  Cape 
Folger  core  taken  from  data  obtained  by  Wakahama 
are  given  in  Budd  (1972).  The  usual  presentation  for 
ice  crystal  orientation  data  is  the  fabric  diagram  which 
gives  the  projection  of  the  cumulative  crystal  orienta- 
tion from  a hemisphere  onto  a plane.  An  example  of 
a fabric  diagram  from  Cape  Folger  was  given  in  Figure 
2 and  the  derived  vertical  distribution  in  Figure  3.  The 
first  moment  of  the  distribution  (mean  angle)  from 
each  depth  is 

ir/2 

8 - V 0n(0)A6.  (26) 

The  variation  of  ? with  depth  is  shown  in  Figure  6a. 
The  graph  can  then  be  used  to  compute  the  dielectric 
constant-depth  variation  from  the  following  approxi- 
mation to  eq  25: 


e(z)  = ex+  0.0037  cos?  (2)  (27) 

where  8 is  the  mean  angle  at  depth  given  in  Figure  6a. 
The  variation  of  Ae  = 0.0037  cos?  with  depth  is  given 
in  Figure  6b. 

As  mentioned  earlier,  this  relationship  does  not 
account  for  the  azimuthal  variation  in  crystal  orienta- 
tion; however,  that  quantity  shows  small  variation  for 
this  particular  core. 

Radio-echo  sounding  data 

There  have  been  three  echo  sounding  traverses  of 
the  Cape  Folger  area  carried  out  by  the  Australian 
National  Antarctic  Research  Expedition  (ANARE):  in 
1967  at  a frequency  of  35  MHz,  and  in  1973  and  1975 
at  a frequency  of  100  MHz.  Unfortunately,  the  1975 
traverse  did  not  take  any  measurements  within  200  m 
of  the  core  site,  the  core  having  been  obtained  in  1969. 
The  bed  topography  appears  to  be  so  rough  in  the  Cape 
Folger  area  that  the  depth  determination  from  the  1975 
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Figure  7.  Dielectric  constant  and  power  reflection  coefficient  vs  depth,  In- 
cluding the  effects  of  density  variations,  bubble  elongation  and  orientation, 
and  crystal  fabric.  The  positions  of  the  internal  layers  observed  by  radio- 
echo sounding  are  shown  (The  amplitude  of  the  observed  echoes  is 

unknown  and  their  position,  as  shown,  on  the  dB  axis  is  arbitrary.) 


traverse  is  100  m deeper  than  the  previous  traverses. 

In  addition,  the  transmitting  pulse  length  used  during 
this  traverse  was  long  enough  to  override  any  return 
echoes  from  depths  less  than  300  m.  Hence,  we  con- 
sider only  the  results  from  the  earlier  work.  The  1967 
traverse  indicated  an  ice  depth  of  345  m and  two 
closely  spaced  layers  at  218  and  235  m.  The  traverse 
in  1973  found  an  ice  depth  of  321  m,  a strong  echo 
at  177  m,  and  a weak  one  at  228  m. 

The  reasons  for  the  differences  between  the  two 
measurements  are  open  to  discussion.  One  possibility 
is  the  different  frequencies  used  in  the  two  years.  If 
the  enhancement  of  a return  is  due  to  an  effect  with 
an  orientational  character,  such  as  bubble  elongation, 
or  crystal  orientation,  then  the  strength  of  the  return 
would  also  depend  on  the  relative  orientations  of  the 
transmitting  and  receiving  antennas.  We  will  take  the 
above  work  to  indicate  the  existence  of  two  layers  at 
Cape  Folger:  one  at  1 77  ± 1 5 m and  a second  one  in 
the  218-  to  235-m  region,  given  by  228  ± 15  m.  The 
ice  depth  from  the  core  closely  corresponds  to  the  ice 
depth  determined  from  the  1 973  traverse,  indicating 
proximity  of  the  locations.  The  errors  on  the  layer 
positions  are  based  on  the  uncertainty  in  reading  the 
travel  times  from  the  A-scan  photographs  of  about 
1 0.2  psec,  leading  to  a possible  layer  depth  error  of 
± 15  m. 

There  might  be  some  objection  to  using  results  from 
the  above  traverses  in  1967  and  1973  because  they 
were  done  at  different  frequencies.  The  frequency 


difference  appears  in  eq  5 in  the  tanhyjrfj  term.  How- 
ever, with  the  choice  we  make  for  d{  = 10  m,  and 
using  Cj  «•  3.17  in  eq  3 foryj,  tanhyj</j  » 1.0,  for 
either  35  MHz  or  100  MHz.  In  this  range  of  frequencies 
and  for  our  choice  of  dit  our  results  should  therefore 
be  frequency  independent. 


RESULTS  AND  CONCLUSIONS 

As  indicated  previously,  the  only  effects  felt  to  be 
of  importance  in  explaining  the  layering  at  Cape  Folger 
are  density  fluctuations,  bubble  elongation  and  bubble 
orientation,  and  fluctuations  in  ice  crystal  fabric  orienta- 
tion. Figure  7 shows  the  dielectric  constant  and  the 
power  reflection  coefficient  as  a function  of  depth 
including  all  these  effects.  ea  is  the  dielectric  constant 
with  the  incident  E field  aligned  along  the  long  axis, 
a,  of  the  bubbles,  eb  the  dielectric  constant  for  the  E 
field  along  the  short  axis,  b,  and  the  dielectric 
constant  for  random  bubbles.  The  case  of  purely 
spherical  bubbles,  eJ(  which  would  show  just  the  effect 
of  density  fluctuations,  has  been  excluded  from  Figure 
7,  since  as  was  noted  previously,  it  does  not  differ 
appreciably  from  e^.  The  location  of  the  internal 
layering  from  the  radio-echo  records  is  also  indicated 
in  Figure  7. 

The  layering  is  seen  to  be  located  very  close  to  the 
position  of  the  highest  values  of  the  reflection  coef- 
ficient, which  would  be  the  predicted  depth  for  enhanced 
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Figure  8a.  The  quantity  elb  and  the  corresponding  reflection  coefficient  vs 
depth  with  and  without  crystal  effects. 
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Figure  8b.  The  quantity  eb  and  corresponding  power  reflection  coefficient 
vs  depth  as  in  Figure  8a. 


reflections.  The  enhanced  reflections  do  appear  to 
be  at  a slightly  lower  depth  than  the  measured  layer- 
ing, but  the  data  points  used  to  calculate  the  permit- 
tivity were  taken  every  10  m,  giving  rise  to  a ± 5-m 
uncertainty  in  our  calculations  of  the  permittivities. 
Given  the  uncertainty  in  the  echo  sounding,  that  the 
internal  layering  differs  from  the  region  of  highest 
reflection  coefficient  by  only  10  or  15  m is  well  within 
the  errors. 

The  dielectric  constants  in  Figure  7 all  display  the 
same  variation  as  a function  of  depth,  though  eb  has 
the  highest  values  of  the  reflection  coefficient  and 
indicates  that  it  would  yield  the  most  enhanced  return. 


However,  the  gross  structure  of  the  dielectric  constant 
is  really  determined  by  the  density  fluctuations,  as  can 
be  noted  by  recalling  that  eib  is  a measure  of  the  density 
effects  alone  and  that  it  indicates  enhanced  reflections 
at  the  same  locations  as  the  other  dielectric  constants. 
While  the  overall  structure  of  the  reflection  coefficient 
is  determined  by  the  density  fluctuations,  the  bubble 
orientation  could  change  the  coefficient  by  a few  dB’s, 
with  eb  providing  an  enhancement  and  e,  a decrease 
over  the  values  for  e^. 

Figure  7 indicates  the  bubble  elongation  and  bubble 
orientation  effects  by  noting  the  differences  between 
ea,  eb  and  e^.  Figures  8a  and  8b  show  the  crystal  fabric 
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Figure  9a.  The  dielectric  constant  and  reflection  coefficient  for 
the  various  cases  of  bubble  configurations  as  a function  of  depth 
using  the  linear  regression  density. 


Figure  9b.  The  dielectric  constant  and  reflection 
coefficient  for  crystal  fabric  effects  as  a function 
of  depth  using  the  linear  regression  density. 


orientation  effects.  In  these  figures  e^  and  eb  are 
plotted  as  a function  of  depth  with  and  without  the 
crystal  effects;  a graph  of  et  would  indicate  the  same 
sorts  of  differences  and  so  is  not  included.  Both  of  the 
figures  indicate  that  changes  in  the  crystal  fabrics  can 
lead  to  a I - or  2-dB  change  over  the  results  without  the 
crystal  information.  The  inclusion  of  the  crystal  effects 
lowers  the  reflection  coefficients  because,  in  eq  27, 
e^_is  the  dielectric  constant  without  crystal  effects,  and 
any  deviation  of  the  c-axis  from  the  horizontal  increases 
the  dielectric  constant  and  decreases  the  reflection  co- 
efficient. 

Figures  9a  and  9b  indicate  further  that  the  structure 
of  the  dielectric  constant  is  determined  by  the  density 
fluctuations.  They  are  a result  of  carrying  out  a linear 
regression  of  the  density  upon  depth.  This  results  in 
the  fit 

pL(/)  = 907.1 +0.0144/  (28) 

where  pL  (/)  is  the  density  in  kg/m3  and  / is  the  depth 
in  meters.  This  line  is  also  plotted  in  Figure  4.  From 
equation  21 , this  density  may  be  used  to  define  a 
"linear  density"  dielectric  constant  for  spherical  bubbles 

eL  = 3.1396  + 0.000044/.  (29) 

This  is  plotted  on  Figures  9a  and  9b;  the  reflection  co- 
efficient for  eL  is  - 89  dB,  irrespective  of  depth  density 
based  on  a "step"  change  in  the  dielectric  constant  at 
10-m  intervals  given  by  (|+(  - e(  = 0.00044.  The  other 


quantities  in  the  figures  are  the  permittivities  and  re- 
flection coefficients  with  the  inclusion  of  the  air  bubble 
effects  (Fig.  9a)  and  crystal  fabric  effects  (Fig.  9b). 
Without  the  density  fluctuations,  all  the  reflection 
coefficients  are  less  than  about  - 79  dB  and  average 
about  - 90  dB.  Thus,  removing  the  density  fluctuations 
removes  the  wide  variation  of  dielectric  constants,  as 
shown  in  Figure  7,  and  also  decreases  the  range  of  high 
reflection  coefficient  from  - 55  to  - 60  dB  to  - 79  to 
- 90  dB  in  the  region  of  the  layering. 

We  note,  however,  that  the  calculated  high  reflection 
coefficient  associated  with  the  layer  at  177  ± 15  m is 
based  on  a single  density  measurement  at  195  m depth. 

If  this  measurement  is  incorrect,  we  see  from  Figure  4a 
that  a relative  peak  in  the  reflection  coefficient  would 
still  exist  based  on  the  air  bubble  orientation  changes 
even  when  the  density  variation  is  smoothed  out  by  a 
linear  regression  line.  Although  the  calculated  reflection 
coefficient  is  some  20  dB  down  from  the  possible  den- 
sity effect,  it  is  still  above  - 80  dB  and  represents  a 
"viable"  reflection  cause.  Its  calculated  position  as 
shown  in  F igure  9a  also  more  nearly  coincides  with  the 
measured  position  of  the  radio-echo  layer,  although 
given  the  errors  on  both  the  calculation  and  measure- 
ment, this  may  not  be  a significant  feature. 

The  bubble  effects  give  a wider  range  than  the  crystal 
effects,  leading  to  the  conclusion  that  the  bubble  effects 
would  then  be  a more  likely  source  for  any  enhance- 
ment than  fabric  changes,  and,  as  shown  in  Figures  8a 
and  8b,  the  crystal  effects  actually  decrease  the  reflection 
coefficient,  while  the  bubbles  aligned  with  theii  short 


II 


axes  along  the  incident  E field  provide  an  enhancement, 
as  shown  in  f igure  7.  The  various  physical  properties 
of  the  Cape  Folger  core  indicate  that  the  internal  layer- 
ing is  most  likely  determined  by  the  ice  density  fluctua- 
tions. Changes  in  bubble  elongation  and  orientation 
provide  the  next  most  important  effect,  while  changes 
in  the  crystal  fabric  are  the  least  important.  Reflections 
in  this  core  caused  by  impurities  are  considered  un- 
likely. 

The  glaciological  implications  of  the  data  are  en- 
hanced by  one  additional  piece  of  evidence.  Recent 
radio-echo  sounding  in  this  area  indicated  that  at  least 
one  internal  reflection  had  a strength  dependent  on  the 
orientation  of  the  receiving  antenna  to  the  surface 
(Budd  and  Harrison,  personal  communication). 

Since  the  radio-echo  strength  of  internal  layering  is 
dominated  bv  the  density  differences,  features  usually 
displaying  an  isotropic  or  non-orientational  character, 
we  must  examine  the  available  evidence  to  see  if  there 
is  an  explanation  consistent  with  an  echo  caused  by  a 
change  in  density  but  also  exhibiting  an  orientational 
or  non-isotropic  dependence.  From  Figures  4 and  5, 
it  appears  that  the  depth  locations  of  the  density  vari- 
ations correspond  well  with  the  most  asymmetric  air 
bubbles,  an  indicator  of  enhanced  deformation  at 
these  depths.  We  conclude  that  the  density  variation 
and  reflection  strength  therefore  reflect  a non-isotropic 
deformational  horizon  at  this  location  rather  than  the 
isotropic  depositional  horizon  usually  associated  with 
density  variations  alone.  The  deformational  event  may 
have  been  a surge  or  strong  localized  shear  layer  be- 
tween rapidly  moving  ice  and  more  stagnant  ice. 

Tracking  of  this  layering,  using  a deformational  assump- 
tion for  occurrence,  may  give  some  insight  into  the 
previous  and  ongoing  internal  mechanics  of  the  ice 
sheet  in  this  region.  We  also  infer  that  the  density  vari- 
ations observed  at  other  locations  in  the  ice  sheet  must 
be  correlated  well  with  deformational  indicators  such 
as  air  bubbles,  fabrics,  and  possibly  grain  sizes  before 
heavy  reliance  is  placed  on  their  use  as  depositional 
markers  in  interpreting  the  prior  history  of  the  ice 
sheets. 
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